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and extraction of the pion form factor at threshold. 
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Abstract 

> . 

Effects due to non-pointlike behaviour of pions in the process e + e" — ► 7r + 7r~7 can 
. arise for hard photons in the final state. By means of a Monte Carlo event generator, 

which also includes the contribution of the direct decay <j> — > 7r + 7r~7, we estimate these 
effects in the framework of Resonance Perturbation Theory. We consider angular cuts 
■ used in the KLOE analysis of the pion form factor at threshold. A method to reveal the 

effects of non-pointlike behaviour of pions in a model-independent way is proposed. 

a: 

& 1 ; 1 Introduction 

& ; 

Final state radiation (FSR) is an irreducible background to the measurement of the hadronic 
cross section with initial state radiation (ISR) events PUI2|- Differently from ISR, whose accu- 
racy is limited by the numerical precision on the evaluation of high order QED diagrams (see, 
for example, [SI Hj and discussion there), the FSR evaluation relies on specific models for the 
coupling of hadrons to photons. Usually the FSR amplitude in the process e + e~ — > 7r + 7r _ 7 is 
evaluated in scalar QED (sQED), where the pions are treated as point-like particles and the 
total FSR amplitude is multiplied by the pion form factor computed in the VMD model [HE]. 
While this assumption is generally valid for relatively soft photons, it can fail for low values of 
the invariant mass of the hadronic system, i. e. when the intermediate hadrons are far off shell. 
In this case possible extensions for FSR, beyond sQED, can be considered. 

As we will show in the following, the most general form for FSR consistent with gauge 
invariance, charge conjugation symmetry and photon crossing symmetry can be expressed in 
term of three form factors, ft, each depending on three independent variables [3[Zj. While 
such a decomposition is general, fi are model dependent. In the paper [Hj the prediction for 
fi in the framework of the Resonance Perturbation Theory (RPT) was considered. RPT is a 
model based on Chiral Perturbation Theory (%PT) with the explicit inclusion of the vector and 
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axial-vector mesons, p (770) and ai(1260) jS|. Whereas gives correct predictions on the 
pion form factor at very low energy, RPT is the appropriate framework to describe the pion 
form factor, and to satisfy QCD high energy behaviour, at intermediate energies (E ~ m p ) jH] 1 . 

We point out that the FSR process, besides being of interest as an important background 
source, is of interest by itself, since it allows to get information about pion-photon interaction 
at low energies. 

In this paper we will present the results obtained by implementing the RPT amplitude for 
FSR into the Monte Carlo event generator for the process e + e~ — > it + -n~-y [T]. We also included 
the rare decay <ft 7r + 7r~7 [9 j in our generator. 

Since most of the effects arising in the FSR are model-dependent, we conclude this paper 
by suggesting a way, based on Monte Carlo, to test possible effects beyond sQED, in a model- 
independent way. 



2 General parametrization for n + 7v 7 final state 

The cross section of the reaction 

e + (pi) + e~(p 2 ) -> rt + (p + )ir~(p^(k) 
with the photon emitted in the final state, can be written as 

da = —1— [ 5\P - P+ -p--k) fP±fP^ \ M \\ (1) 
2s(2vr) 5 J V ; 8E + E_u 1 1 ' V ; 

where P = p\ + p 2 , s = P 2 and 

M = ~M^(- Pl ) 7/ ^(p 2 )e*. (2) 

The tensor Mp 1 describing the process J*(P) — > K + (p + )-K~{pJ)^{k) is model-dependent. 
However, based on charge-conjugation symmetry, photon crossing symmetry and gauge invari- 
ance it can be expressed by three gauge invariant tensors (see Appendix A in Ref. [6]): 

M^(P, k, I) = -ie\r^h + T?h + rffs) = -ie 2 M^{P, k,l), l=p+-p- t (3) 
7f = k"P v -g^k- P, 

T f = k ■ l(FP u - gTk ■ I) + l v {k»k -l-Fk- P), 
T £ v = p 2 ( g ^k ■ I - m v ) + P»(l v k -P-P v k-l). 

While the decomposition (jSJ) is general, the exact value of the scalar functions /, (form 
factors), each depending in terms of three independent variables, are determined by the specific 
FSR models. 

In sQED for the functions we have j3] 

rsQED __ 2fc " P psQED _ rsQED __ n , ,s 

Jl ~ (k ■ P) 2 - (k ■ I) 2 ' 12 ~ (k ■ P) 2 - (k ■ I) 2 ' 73 " ' K ' 



1 In that paper it was shown that the coupling constants of the effective chiral lagrangian at the order p 4 are 
essentially saturated by meson resonance exchange. 
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Because of Low's theorem, these equations imply that for k — > we have 

2fc ■ PP^P 2 ) -2P,(P 2 ) . n 

hm^o/i = (fc . p)2 _ (fc . 02 . I™*-*/* = (fc . P) 2_ (fc . 2 » hm ^o/s = 0, (5) 

where P^ is the pion form factor describing the interaction 7* — > 7r + 7r~. Thus, for soft photon 
radiation, the FSR tensor is expressed in term of one form factor F n (P 2 ), but in general we 
have three independent form factors describing the FSR process. 
It is convenient to rewrite the form factors fa as 

fi = ft QED + A/,, (6) 

where the functions A/j are the contributions to the form factors beyond sQED, and are 
determined by FSR model. 

2.1 RPT contribution to FSR 

The functions A/j have been calculated in the framework of RPT in 



P 



/ 2 V m p m 2 - P 2 - im p r p (P 2 ) 

(fc • /) 2 (P 2 + fc • P)[4m 2 - (P 2 + I 2 + 2k ■ P)] 
+ D(l)D(-l) + 8D(l)D(-l) 



(7) 



P 2 4m 2 -(P 2 + Z 2 + 2fc-P) 
72 8D(l)D(-l) ' 1 J 

A/3 = D(0=m:-(g» + ? + 2*g + 4«)/4 (9) 

For notations and details of the calculation we refer a reader to jH]. F v , Gy and Fa are 
parameters of the model. The parameters Fy, Gy as well as m p have been estimated by a fit 
to the pion form factor from e + e~ — > 7r + 7r~ data (see next Section). For ai meson we take 
m a = 1.23 GeV and Pa = 0.122 GeV that corresponds to the mean value of the experimental 
decay width V{a x -> vr 7 ) = 640 ± 246 keV jTU]. 



2.2 Pion form factor in RPT 

The pion form factor, that describes p — u mixing and includes also the first excited p'-meson 
state, can be written as: 

FAi) = 1 + (1 - ^Mf)) + ^f^te 2 ), (10) 



where 



r ^ m 2 — g 2 — im r r r (g 2 ) ' ^ 

g 2 is the virtuality of the photon, f n = 92.4 MeV and the value U puJ describes p-u mixing (see 
below). More detailed description of the pion form factor can be found elsewhere |TT| I12j. 
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An energy dependent width is considered for the p and p' mesons: 

3/2 




g2 .~ 4m -e(a 2 



W) =r M/^l t ; 3 _^ 1 -0(^-4^), (12) 



while for the cu-meson a constant width is used, = 8.68 MeV, and m w = 782.7 MeV. We 
assume the parameter IT^, that determines p-u mixing, is a constant and we relate it to the 
branching fraction Br(cu — > 7r + 7r~): 

Br(u 7T + 7T-) = (13) 
1 pi ujTTLp 

The fit of our parametrization to the pion form factor from the CMD-2 data gives |12| : 

m p = 774.97 ± 1.4 MeV, n pw = -2774 ± 291.2 MeV 2 , 
T p = 145.21 ± 2.6 MeV, F v = 154.22 ± 0.5 MeV 

and 

m' p = 1.2 ± 0.2 GeV, r p / = 400 ± 160 MeV, Fyi = 13.19 ± 18.59 MeV, 

x */u = 0.853. Then G v = 64.6 ± 0.3 MeV and Br(u -> vr+vr-) = (0.96 ± 0.19)%. Expanding 
the pion form factor in the region s < 0.35 GeV 2 as in [Tfj : 

F w (g 2 ) ~l+ Pl -g 2 +p 2 -g 4 , (14) 

gives the following values: p 1 = 1.15 ± 0.06 GeV" 2 , p 2 = 9.06 ± 0.25 GeV" 4 , y?jv ~ 0.13. 



2.3 Scalar contribution to FSR 

At s — rr& an additional contribution to the final state 7r + 7r~7 is given by the direct rare decay 
cf) — > 7T + 7r _ 7. As it was shown in [HIE] this process affects the form factor fi of Eq.(j3J): 

if? 

Mct> = —F^Q 2 )u(- Pl )^u(p 2 )e:rr, (15) 



The direct decay is assumed to proceed through the / intermediate state: <fi — > / 7 — *■ 
7r+7r - 7, and its mechanism is described by a single form factor f^,(Q 2 ). 

We consider the model described in [H)] where the — > / 7 decay amplitude is generated 
dynamically through the loop of charged kaons. The form factor reads: 

f K+K- {n 2\ = 9<j,K+K- gf n+n- 9f K+K- j ( Q^_\ iS B (Q 2 ) /-| fi \ 

J * W> 2^mUml-Q 2 + ReU f0 (ml)-U f0 (Q^) Wk^k) ' 1 ' 



where /(.,.) is a function known in analytic form [5J El and Sb(Q 2 ) = bWQ 2 — 4m 2 ., b = 
75°/GeV [THj . The term i?eil/ (m 2 ) — II/ (Q 2 ) takes into account the finite width corrections 
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Figure 1: Contribution to the FSR cross section dopjdQ 2 in the region 0° < # 7 < 180°, 
0° < 0-k < 180°. RPT is represented by circles, sQED by crosses, <fi by triangles, while the 
difference between RPT and sQED is indicated by squares. (Left) corresponds to s = m|, 
(right) to s = 1 GeV 2 (i.e. below the resonance), where the <fi resonant contribution has been 
amplified by a factor 100. In (left) the results from our event generator are compared with the 
analytic calculation, shown by solid line. 

to the fo propagator ^B]. A fit to the KLOE data <fi — > 7r°7r°7 2 gives the following values of 
the parameters jT9"] : 

m fo = 0.962 GeV, 9foK+K ~ = 1.29 GeV 2 , 9f f +K ~ = 3.22. (17) 

4:% qi . _ 

A refined version of this model includes the presence of the a meson in the intermediate 
state [IH1 EH EHj • Such an extension of the model improves the description of the data at low 
Q 2 and will be considered in a forthcoming paper ^2]- Therefore, for the following results, 
we will only consider the presence of the fo meson. This will not affect our proposal, to be 
discussed in Sec. 4. 

2.4 Other contributions 

We included in our program the channel 7* — > p ± 7r =F — > 7r + 7r _ 7, whose amplitude has been 
evaluated in RPT model. However, in agreement with the calculation given in [6J, we found a 
negligible contribution of this channel and, for the sake of simplicity, we discard the effects on the 
following results. We also did not include the contribution due to two-resonance intermediate 
states as, for example, — > pn, which was found to be negligible in the 7r°7r°7 final state [TH] . 

2 r(/ - 7T+7T-) = |r(/ - TTTT) 
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Figure 2: The ratio dax/dai as function of the invariant mass of the two pions, in the region 
50° < # 7 < 130°, 50° < 9 n < 130°, for different models of FSR. (Left) refers to s = mj, (right) 
to s = 1 GeV 2 . 



3 Numerical results 

The differential cross section for the reaction e + e~ — > 7r + 7r~7, where the FSR amplitude (Mfsr) 
receives contributions both from RPT (Mrpt) and the — > 7r + 7r _ 7 decay (M^) can be written 
as: 

(Lot ~ + Mfsr\ 2 = ^cr/ + d<JF + dojF-, (18) 

do/ ~ \M ISR \ 2 , 

da F ~ |iVW| 2 + |M/ + 2Re{ilW-M;}, 
da JjF ~ 2Re{M /5i? • (ilW + M^)*}. 



The interference term daip is equal to zero for symmetric cuts on the polar angle of the pions pQ. 

The different contributions to the FSR differential cross section dap, evaluated at s = m^, 
are shown for in Fig.l, left, for the full angular range 0° < # 7 < 180°, 0° < 9 n < 180°. A 
good agremeent between the results of the Monte Carlo simulation (points), with the analytic 
prediction (solid line) is found. It can be noted that at low Q 2 the resonant contribution 
(i.e. the term proportional to | Af^| 2 in Eq. (JUJ)) is quite large and, therefore, the additional 
contribution beyond sQED, can be revealed only in the case of destructive interference between 
the two amplitudes (Re(M^pr ■ Mi) < 0). Published data from KLOE experiment [2] are in 
favour of this assumption, which we will use in the following. 

In Fig. 2 we show the values of dar/dai for the angular cuts of the KLOE large angle analysis 
50° < # 7 < 130°, 50° <9 W < 130° [21 , with and without contributions from RPT and direct 
decay, for a hard photon radiation with energies E 1 > 20 MeV. 

Three distinctive features can be noted: (1) the peak at about 1 GeV 2 corresponds to the 
fo intermediate state for the — > tttt^ amplitude; (2) the presence of RPT terms in the FSR 
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Figure 3: Ratio of FSR cross section in the framework of RPT, respect to sQED, when the < 
contribution is (or not) taken account. The angular region is 50° < 6 1 < 130°, 50° < 9 W < 130 c 
(Left) is for s = rr&, (right) is for s = 1 GeV 2 . 



are relevant at low Q 2 , where they give an additional contribution up to 40% on the ratio 
d&RPT+<t>/ ' da S QED+4>i (as shown in Fig. 3, left); (3) the negative interference with the direct 
decay amplitude reduces dap and its dependence on FSR model at low Q 2 (see Fig. 2, left, 
down, and Fig. 3, left). 

Fig. 4, left, shows the effects of RPT and terms in the forward-backward asymmetry. As 
expected the presence of RPT gives relevant effects at low Q 2 region, while the presence of a 
bump at high Q 2 is due to the direct decay. 

In order to reduce the background from — > 3n decay on the measurement of the pion form 
factor at threshold, KLOE has taken more than 200 pb _1 of data at 1 GeV [221 • In this case, 
as shown in Fig. 1, right, the resonant contribution is suppressed (dap with and without the 
direct decay almost coincide), see also Figs. 2 and 3, right. Therefore the main contribution 
beyond sQED to the FSR cross section and the asymmetry (see Fig. 4, right) comes from RPT. 

4 Model-independent test of FSR and extraction of pion 
form factor at threshold 

Contributions to FSR beyond sQED, as in the case of RPT, can lead to sizeable effects on the 
cross section and asymmetry at threshold, as shown in Figs. 3 and 4. Precise measurement 
of the pion form factor in this region needs to control them at the required level of accuracy. 
This looks like a rather difficult task, if one thinks that effects beyond sQED, as well as the 
contribution from — > 7r + 7r~7, are model dependent. 

One can think to construct a general amplitude for the e + e~ — > 7r + 7r~7, according to some 
underlying theory, and try to determine the free parameters by a constrained fit on specific 
variables (like mass spectrum, charge and forward-backward asymmetry, angular distribution, 
etc.). Particularly for the charge asymmetry, it has been proved to be a powerful tool to 
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Figure 4: Forward-backward asymmetry, in the kinematical region 50° < # 7 < 130°, 50° < 6 W < 
130°, when RPT and <fi contributions are included, compared with the asymmetry calculated 
in sQED. (Left) is for s = m 2 , (right) is for s = 1 GeV 2 . 

discriminate between different models of — > tt + tt~^ However when the number of the 
parameters is large, correlations between the parameters of the model can arise and spoil the 
effective power of these fits. The situation becames even worse if also the pion form factor has 
to be extracted with the same data. As an example, in the case of RPT model, if we consider 
only p and u contribution to the pion form factor and the p and a± contribution to FSR the 
number of free parameters is already six. The presence of the </> direct decay adds additional 
free parameters. 

The possibility to determine some of the parameters by external data can strongly help, 
as in the case of the — > 7r + 7r~7 amplitude, which can be determined by the 7r°7r°7 channel 
copiously produced at DA$NE. An additional source of information which will be used to 
determine the contributions to FSR beyond sQED in a model-dependent way, is the energy 
dependence of the FSR amplitude on the e + e~ invariant mass s. 

Let us write the differential cross section for the emission of one photon in the process 
e + e~ — > 7r + 7r~7 as function of the invariant mass of the two pions: 

daT ^ -'^(Q 2 )i 2 ^(g 2 ) + fe) , (19) 

where H S (Q 2 ) is the so called radiation function, which accounts for ISR emission, and (j^r^J 

is the differential cross section for the emission of a photon in the final state. We explicity 
put in evidence the dependence of each quantity on the e + e~ invariant mass (s). Since we will 
consider only symmetric angular cuts for pions, the interference term between initial and final 
state radiation has been neglected. 

The FSR differential cross section, {^§^j , is dominated at relatively high Q 2 by the 



dQ 2 Js 1 " ^ ' \dQ 



S 



contribution coming from sQED (M s q ED ) and direct decay (M^): 

(^fP) s HM sQSD + M/. (20) 

Contributions beyond sQED (AM) are expected to be important at low Q 2 . They introduce 
an additional term (AM) in the above expression: 

i^)s ~ \M sQED + AM + M (p \ 2 = (21) 

= \M sQED + M^\ 2 + \AM\ 2 + 2Re{AM ■ (M sQED + M^)*}. (22) 
We will now consider the following quantity: 

/ dcr T \ / dfjQ_ED+±\ 

Y S (Q 2 ) = h = \ F *(Q 2 )\ 2 + AF *(Q 2 )> ( 23 ) 

where AF S ~ (|AM| 2 + 2Re|AM • (M sQED + Afy)*}) /#,. 

If no contribution beyond sQED is present (AM = 0), Y S (Q 2 ) coincides with the square of 
the pion form factor, independently of the energy at which it is evaluated, while any 
dependence on s is only due to additional contribution to FSR beyond sQED. In particular, 
the difference of Y S (Q 2 ) computed at two beam energies (si and s 2 ), can only come from FSR 
beyond sQED: 

AY{Q 2 ) = Y S1 (Q 2 ) - Y S2 (Q 2 ) = AF Sl (Q 2 ) - AF S2 (Q 2 ) (24) 

Therefore, before extracting the pion form factor at threshold, we suggest to look at the differ- 
ence AY(Q 2 ), which can be used to estimate the contribution beyond sQED to FSR amplitude 
in a model independent way. 

As realistic application of this procedure, we consider the case of DA$NE, where KLOE 
has already collected more than 200 pb -1 at 1 GeV 2 and 2.5 fb -1 at m 2 , which, in the range 
Q 2 < 0.35 GeV 2 , correspond to O(10 3 ) and O(10 4 ) events respectively in the region 50° < 
# 7 < 130°, 50° <9 W < 130°. We will consider RPT as model for the effects beyond sQED. 

Fig. 5, left, shows the quantity Y S (Q 2 ) at s\ = 1 GeV 2 and at s 2 = m 2 ^ when no additional 
RPT term is included in FSR. As expected, each of these quantities coincides with the square 
of the pion form factor \F n (Q 2 )\ 2 , shown by solid line. The difference AY(Q 2 ) is shown in 
Fig. 5, right, which is consistent with zero as expected. A combined fit of Y S (Q 2 ) to the pion 
form factor (see Eq. (JHJ)) gives the following values: p\ = 1.4 ± 0.186 GeV -2 , p 2 = 8.8 ± 0.73 
GeV -4 , x 2 / u = 0.25, in agreement with our results at the end of Sec. 2.2. 

A different situation appears if FSR emission from pions is modeled by RPT. In this case, as 
shown in Fig. 6, right, the difference AY(Q 2 ) ^ and the quantities Y S (Q 2 ) cannot be anymore 
identified with \F n (Q 2 )\ 2 , (see Fig. 6, left) 3 . Prior to the fit of the pion form factor, in this case, 
such additional contribution must be understood. 

Before concluding, we would like to point out the main advantages of our proposal: 

• NLO correction to ISR (as multi-photon emission) can be computed by Monte Carlo and 
included in H s ; 



3 Destructive interference between RPT and (j> — > tt + it 7 amplitudes tends to cancel out the effects beyond 
sQED at s = m 2 , (see Fig. 2, left). Therefore the quantity Y S (Q 2 ) almost coincides with the pion form factor. 
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Figure 5: Left: Y S (Q 2 ) at s = 1 GeV 2 (triangles), and at s = rr& (circles), when FSR includes 
only sQED and contribution. The pion form factor \F n (Q 2 )\ 2 is shown by solid line. Right: 
The difference AY(Q 2 ). 




Figure 6: Left: Y S (Q 2 ) at s = 1 GeV 2 (triangles), and at s — m 2 ^ (circles), when FSR includes 
RPT and (f) contribution. The pion form factor \F n (Q 2 )\ 2 is shown by solid line. Right: The 
difference AY(Q 2 ). 
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• The quantity y— ^r^J is an input parameter of our procedure, and can be computed 
numerically by Monte Carlo; 

• The amplitude for — > 7r + 7r~7 is taken from the 7r°7r°7 channel, therefore its description 
is not restricted to the presence of /o only, as in our simulation; 

• A clear advantage of procedure based on a Monte Carlo event generator is that it allows 
to keep control over efficiencies and resolutions of the detector and fine tuning of the 
parameters. 

Even if the main limitation of the method could come by the uncertainty on the parameters 
of — > 7r + 7r~7 amplitude, expecially at low Q 2 , we believe that the new data on — > 7r°7r°7 from 
KLOE will allow a precise description of this amplitude. In any case, in agreement with [20] 
we strongly recommend to check the amplitude by using charge asymmetry and to compare 
with spectrum of the 7r + 7r~7, at least at high Q 2 , where the pointlike approximation is safe (as 
done in j2~3]). 

5 Conclusion 

Test of FSR at threshold in the process e + e~ — > 7r + 7r~7 is a rather important issue, not only 
for the role of FSR as background to the measurement of the pion form factor, but also to get 
information about pion-photon interaction when the intermediate hadrons are far off shell. At 
s = m| an additional complication arises: the presence of the direct decay <fi — > 7r + 7r~7 whose 
amplitude and relative phase can be described according to some model. By means of a Monte 
Carlo event generator, which also includes the contribution of the direct decay <p 7r + 7r _ 7, we 
estimate the effects beyond sQED in the framework of Resonance Perturbation Theory (RPT) 
for angular cuts used in the KLOE analysis of the pion form factor at threshold. We show that 
the low Q 2 region is sensitive to the inclusion of additional terms in the FSR amplitude given 
by the RPT model. We propose a method which allows to estimate the effects beyond sQED 
in a model-independent way. We found that the deviation from sQED predicted by RPT can 
be observed with the current KLOE statistics. 
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